A Low CMB Quadrupole from Dark Energy Isocurvature Perturbations 
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We explicate the origin of the temperature quadrupole in the adiabatic dark energy model and 
explore the mechanism by which scale invariant isocurvature dark energy perturbations can lead to 
its sharp suppression. The model requires anticorrelated curvature and isocurvature fluctuations and 
is favored by the WMAP data at about the 95% confidence level in a flat scale invariant model. In 
an inflationary context, the anticorrelation may be established if the curvature fluctuations originate 
from a variable decay rate of the inflaton; such models however tend to overpredict gravitational 
waves. This isocurvature model can in the future be distinguished from alternatives involving a 
reduction in large scale power or modifications to the sound speed of the dark energy through the 
polarization and its cross correlation with the temperature. The isocurvature model retains the 
same polarization fluctuations as its adiabatic counterpart but reduces the correlated temperature 
fluctuations. We present a pedagogical discussion of dark energy fluctuations in a quintessence and 
k-essence context in the Appendix. 



I. INTRODUCTION 

The first data release of the WMAP data [lj was on 
the whole in spectacular agreement with the ACDM 
model 0, 13- However, as originally discovered by 
COBE, the measured value of the quadrupole seems 
low compared to the model prediction. Originally it 
was estimated that the probability of measuring such 
a low or lower quadrupole was 0.7% assuming that the 
ACDM model was correct 0. However, subsequent stud- 
ies with more detailed methods of dealing with fore- 
grounds have estimated this probability to be closer to 
4% [H Ga, p. H, lal . The alignment of the quadrupole and 
octopole JM l9l llOl IX ll and various asymmetries in the 
data [l3. fla. fl4L fl5l fla. Il7| have also being considered. 

This probability of such a low or lower quadrupole 
would not be particularly anomalous if the quadrupole 
were treated on par with all other multipoles. Such points 
would be expected to occur just by chance in such a large 
data set. In fact there are several other equally or more 
anomalous multipoles. However the quadrupole is par- 
ticularly intriguing in that it represents the largest ob- 
servable angular scale and so may be a good probe of new 
physical effects. In particular, it is the multipole whose 
power has the most significant contribution from length 
scales that are on and even above the horizon at dark 
energy domination. 

Various explanations for the low quadrupole have been 
proposed. They include a cut off in the primordial power 
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l31l |. a small universe |2l l32| |33, 34, 35] and perturbations 
in the dark energy J33, H3, OH, H^, E3 • The cut-off and 
small Universe models work by reducing the Sachs Wolfe 
effect in the quadrupole. The dark energy perturbation 
models work by modifying the Integrated Sachs Wolfe 
effect from the dark energy. In this paper we focus on 
the latter class and in particular a generalization of the 
correlated isocurvature model introduced by |3q . 

We begin in <|n]with a general discussion of the origin 



of the temperature quadrupole in the adiabatic model 
and explore its relationship to the low multipole polariza- 
tion. In illlll we show how the properties of the adiabatic 
quadrupole point to a specific class of isocurvature mod- 
els that can cancel the Sachs- Wolfe contributions to the 
quadrupole. We compare and contrast this model with 
alternate solutions and show that the polarization will 
be in the future a useful discriminator. In fll VI we assess 
the likelihood of substantial isocurvature perturbations 
in light of the WMAP data. We discuss an inflationary 
context for such perturbations in Wl but show that in the 
simplest models gravitational waves are over-predicted. 
We also include a pedagogical Appendix on dark energy 
perturbations in a quintessence and k-essence context. 



II. QUADRUPOLE TRANSFER FUNCTION 

In an adiabatic model with dark energy, the CMB tem- 
perature quadrupole receives its contributions from two 
distinct effects: the (ordinary) Sachs- Wolfe (SW) effect 
from temperature and metric fluctuations near recombi- 
nation and the Integrated Sachs- Wolfe (ISW) effect from 
changes in the metric fluctuations due to the dark en- 
ergy. These effects are quantified by the CMB tempera- 
ture transfer function. 

Let us define the two dimensional CMB transfer func- 
tions as the mapping between the power in the initial 
curvature fluctuations Q in the comoving gauge [see Ap- 
pendix, Eqn. l|A"3"5)l ] 
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where X, X' £ <d,E the temperature fluctuation and E- 
mode polarization respectively. 
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FIG. 1: The quadrupole transfer function in the fiducial adi- 
abatic model (see text). The temperature quadrupole receives 
Sachs- Wolfe (SW) contributions peaking around k = 0.0002 Mpc -1 
and Integrated Sachs-Wolfe (ISW) contributions peaking around 
k = 0.001 Mpc" 1 but extending to k ~ 0.01 Mpc" 1 . The ISW 
effect arises from the dark energy dominated regime z < 1. The 
polarization arises through rescattering of quadrupole anisotropies 
at z > 1 and hence reflects the SW contributions. The cross corre- 
lation is proportional to the product of the transfer functions. 



We employ a comoving gauge Boltzmann hierarchy 
code |41j for numerical solutions of the transfer function. 
These are shown for the temperature and polarization 
quadrupole in Fig. ^ Here we have chosen fiducial values 
for the cosmological parameters that are near the maxi- 
mum likelihood model from WMAP: a dark energy den- 
sity relative to critical of Qq = 0.73, non-relativistic mat- 
ter density Vt m h 2 = 0.14, baryon density fl^h 2 = 0.024, 
optical depth to reionization r = 0.17, dark energy equa- 
tion of state wq = Pq/pq = —1 in a spatially flat 
universe. The resulting temperature power spectrum is 
shown in Fig. [21 compared with the WMAP data for a 
scale invariant spectrum of initial perturbations 
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where S& = 5.07 x 10~ 5 and the tilt n = 1. 

Notice that the temperature contributions from the 
SW and ISW effects are comparable in magnitude and 
well-separated in scale. Hence they add nearly in quadra- 
ture in Eqn. |2"|. In the polarization, the quadrupole 
comes mainly from the SW effect as can be seen from its 
dependence on the redshift of reionization or equivalcntly 
the cumulative contributions from z < z max (see Fig.^>). 
The polarization quadrupole gets nearly no contributions 
from 2 < 1 when the dark energy dominates. These prop- 
erties are the key to understanding how to construct a 
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FIG. 2: Temperature power spectrum in the fiducial mod el com- 
pared with the data from the first year WMAP release 1421 . Bands 
arc 68% and 95% cosmic variance confidence regions (see text). 



model with certain desired properties at low temperature 
and polarization multipoles. 

Since the quadrupoles are dominated by large scale 
or low-fc fluctuations, it is useful to examine the origin 
of these properties with a low-fc approximation to the 
transfer functions. In this limit, the transfer function 
is determined by the Newtonian temperature monopole 
G, gravitational potential \P and curvature fluctuation 
$ [see Eqn. QA41J I. i|A42|) for the correspondence to the 
comoving gauge] 



T?(k,vo) 



VW+I) 



(e* + **)ji(kD,) 



cfy(* - *)ji(kD) 



y/2l(l+l) 



:Cde(kD m ) 



+2/ dr)$j e (kD) 

JVm 



(4) 



where the subscripts denote evaluation at recombination 
for " * ", "0" for the present and "m" for some arbi- 
trary time well after radiation domination but well before 
dark energy domination and overdots represent deriva- 
tives with respect to conformal time rj = J dt/a. Here 
we have assumed a spatially flat universe where the co- 
moving distance D = r/o — rj. Since $ w — $ when the 
anisotropic stress is negligible [see Eqn. 1A42|> ]. we loosely 
refer to either as the gravitational potential. 

The two terms in the second line are the SW and 
ISW effects respectively. This approximation accounts 
for the small evolution in the gravitational potential be- 
tween recombination and full matter domination, some- 
times called the "early" ISW effect. Since D m rs £>*, 
this effect adds coherently with those at recombination. 
In the fiducial cosmology, D* w 14Gpc. Since je(x) peaks 



0.5 



-0.5 

-1 
1 

0.5 




-0.5 



(a) £=0.0002 Mpc- 1 




(b) £=0.001 Mpc 



plies |43 



0.01 0.1 

a 



FIG. 3: Time evolution of gravitational potential <1> and 
quadrupole T® in the fiducial adiabatic model compared an adia- 
batic plus isocurvature model with Si = &Qi/C,i = — 15 (A&I; see 
flTTl . (a) k = 0.0002 Mpc -1 SW dominated mode. In the adia- 
batic model, the decay of the potential at z < 1 has little effect on 
the quadrupole due to the inefficiency in the transfer represented 
by J2- A much larger change of A<I> ~ — Q from the isocurvature 
perturbation can cancel the SW quadrupole. (b) k = 0.001 Mpc -1 
ISW dominated mode. The adiabatic potential decay transfers ef- 
ficiently onto the quadrupole moment. The additional isocurvature 
perturbations decay by dark energy domination and have little ef- 
fect (see Fig.[|}. 



at x « £ + 1/2, the SW contributions peak at x = 5/2 
or k ~ 0.0002 Mpc in Fig. ^ The temporal evolution for 
the potential and quadrupole are shown for this mode in 
Fig. Hi- 

The ISW effect arises from the change in the gravita- 
tional potential due to the dark energy. In the adiabatic 
model, the smooth dark energy density accelerates the 
expansion without enhancing the density perturbations 
and leads to a decay of the potential at z < 1 (see Fig. 121. 
Naively, one would suppose that a decay of A<& = — £i/10 
would be sufficient to cancel out the SW effect of d/5. 
For superhorizon scaled adiabatic fluctuations in a flat 
universe the conservation of the comoving curvature im- 
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where the approximation assumes a m <^C 1 is some epoch 
near the beginning of matter domination so that contri- 
butions near the lower limit of the integral may be ig- 
nored. Here H — a~ l da/dt is the Hubble parameter. In 
the fiducial model A<I> ss — 0.14^. However the distance 
to z < 1 is much smaller than D* and that degrades 
the efficiency with which the ISW effect contributes in 
Eqn. (0J. The efficiency factor J2(kD) implies that one 
requires factor of 10 greater change in the potential (or 
A$ ~ — £;) to affect the quadrupole substantially (see 
Fig. EJ). At the peak of the SW effect in k, the ISW effect 
has little effect (see Fig.0. 

Contributions from the ISW effect in the quadrupole 
actually originate from scales where kD^c ~ 5/2 where 
"de" denotes the epoch of dark energy-matter equality. 
Because the distance changes rapidly with redshift lo- 
cally, the ISW effect is spread out across over a factor of 
ten in physical scale with a peak centered near k = 0.001 
Mpc -1 (see Fig. ^,). At this wavenumber J2(kD) has a 
peak near z = 1 corresponding to an efficient transfer of 
power. Moreover since the decay occurs on the expan- 
sion time scale the oscillations in k from ji{kD) in the 
ISW integrand of Eqn. Q are washed out in the trans- 
fer function. Physically this reflects the cancellation of 
radial modes as photons travel in and out of decaying 
gravitational potentials along the line of sight. 

To lower the predicted value of the quadrupole, one 
can alter the fiducial model to lower the SW effect, the 
ISW effect or both. Since both effects contribute nearly 
equally reducing one or the other can at best halve the 
power. Of course, due to cosmic variance, it is possible 
that the observed quadrupole results from a lack of angu- 
lar power in our given realization of the fiducial model. 
In Fig. |21 we show the 68% and 95% cosmic variance con- 
fidence regions assuming that Cf e is distributed as a \ 2 
with 2^+1 degrees of freedom around the fiducial model. 
However again, a simple lack of power on large physi- 
cal scales for our last scattering (recombination) surface 
is not sufficient. Unless our local volume also lacks in- 
termediate scale power as well, a chance occurrence of 
a low observed quadrupole would result from a chance 
cancellation of the SW and ISW effects. 

The low multipole polarization and cross spectra pro- 
vides key additional information to discriminate between 
alternatives. In the large scale limit, it is approximately 
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where r is the Thomson optical depth as measured from 
the observer. Consider the transfer function of the po- 
larization quadrupole. Fig. ^ shows that in the fiducial 



model it is generated before the temperature quadrupole 
is modified by the ISW effect. Consequently, the transfer 
function also peaks near the large scales of the SW peak 
of k ~ 0.0002 Mpc -1 . As Eqn. © shows this overlap 
is also the origin of the temperature-polarization cross 
correlation. In Fig. we show the EE and QE power 
spectra of the fiducial model. 

If the explanation of the observed low quadrupole in- 
volves the dark energy, either through a dynamical ef- 
fect or chance cancellation, one would expect a polariza- 
tion quadrupole and hence EE power that is not anoma- 
lously low compared with the fiducial model. If on the 
other hand it involves an actual lack of predicted long- 
wavelength power in the model or by chance, both spec- 
tra would be low. Finally, if the explanation involved 
only the reduction of the ISW effect and no modification 
of long wavelength power, then the predicted QE cross 
power spectrum would also remain unchanged. 



III. DARK ENERGY MODELS 

The fact that in the fiducial adiabatic model the tem- 
perature quadrupole receives comparable contributions 
from recombination and dark energy domination through 
the SW and ISW effects raises the possibility that the low 
quadrupole originates in the dark energy sector. In the 
Appendix we present a detailed treatment of perturba- 
tion theory in general dark energy models which provides 
the basis for results in this section. 

The essential element that defines the ISW effect in the 
fiducial model is that the dark energy remains smooth out 
to the horizon scale and hence does not contribute den- 
sity fluctuations to the gravitational potential. In general 
there are two ways to alter this conclusion: modify the 
dynamics of the dark energy so that dark matter fluctua- 
tions remain imprinted on the dark energy or modify the 
initial perturbations in the dark energy sector. Since the 
dark energy has made a negligible contribution to the net 
energy density until recently, the latter represents con- 
tributions from an isocurvature initial condition. That 
dark energy isocurvature conditions can help to lower 
the quadrupole has recently been shown 36]. Here we 
present a general discussion on the requirements of such 
a model. 

The first requirement is that an initial dark energy 
perturbation must survive evolution in the radiation and 
matter dominated epoch and must remain correlated 
with the perturbations in the dark matter. The latter 
condition is required for the dark energy perturbations 
to cancel the adiabatic ones. 

Let us take the dark energy to be a scalar field Q 
with the canonical kinetic term and a potential Vq, i.e. 
quintessence. We treat the more general case of k-essence 
in the Appendix. Given that a quintessence field has an 
effective sound speed c e = 1 [see Eqn. (|A28|) ]. coher- 
ence well inside the horizon and hence cancellation with 
the adiabatic ISW effect is not possible. The dark en- 



ergy isocurvature mechanism then must operate on large 
scales to cancel the SW effect. 

Recall that to achieve a coherent cancellation of the 
SW effect in the quadrupole one requires either a change 
of A<J> ~ — £i/10 in the gravitational potential early on 
when D sa D* or a larger change A<I> ~ — Q at z < 1 to 
compensate for the inefficiency of the transfer of power 
to the quadrupole. Given that observations require that 
wq ~ — 1 today, the former possibility is excluded unless 
wq evolves substantially from its present value. 

In a flat universe the comoving curvature evolves only 
in response to stress fluctuations in the combined or to- 
tal ("T") stress energy tensor of the components (see 
Eqn. |A39|) 
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where the approximation assumes wq sw — 1, a 3> a m , 
and the radiation and hence the anisotropic stress is neg- 
ligible. The generalization of Eqn. JSJ for the evolution 
in the Newtonian potential is [see e.g. [44| Eqn. (52)] 
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Thus it requires a substantial pressure fluctuation to 
make an order unity change to gravitational potential 
during the dark energy dominated regime 
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Note that in the comoving gauge, adiabatic density fluc- 
tuations scale as 6 ~ (fc??) 2 Ci [ see Eqn. (|A46|l ] and are 
negligible outside the horizon. Furthermore we shall see 
below that the order unity coefficient in front of Q is in 
practice substantially greater than unity. 

This requirement severely limits the range of 
quintessence models which can affect the quadrupole. 
As shown in the Appendix, aside from transient initial 
condition effects, an isocurvature perturbation to the 
quintessence field 5Q at best remains constant outside 
the horizon and hence one requires a large initial fluctu- 
ation to the quintessence field. 

A constant superhorizon quintessence field fluctuation 
generically occurs if the background field itself is nearly 
frozen by the Hubble drag so as to only experience a 
range in the potential where 
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can be approximated as constant. More specifically, we 
require that the field not be in the tracking regime or 



(|45j see also Appendix) 
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To see the consequences for the energy density and pres- 
sure, note that aside from a transient decaying mode the 
Klein-Gordon equation l|A33[l has the solution 
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where we have assumed an epoch during which wt = 
Pt/pt the equation of state of the background is con- 
stant. The dark energy density is the sum of the kinetic 
and potential components 
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Since H 2 ex pr decreases with the expansion, if the field 
is potential energy dominated today [wq ~ — 1) then 
it is potential energy dominated for the past expansion 
history. The combination of Eqn. (|12fl and (|13|) shows 
that this will be satisfied if the potential satisfies 



1 

2 

or equivalently with 



11 
Vq 



Vq 



< i. 



e Q 



16ttG \V q 



VA 



and 



fir 



8ttGVq 



the condition becomes 



(14) 



(15) 



(16) 



FIG. 4: Isocurvature dark energy density perturbation evolution 
for a nearly frozen background quintessence field (see text) . On su- 
perhorizon scales the comoving gauge density perturbation remains 
frozen at nearly its initial value. On subhorizon scales, the pres- 
sure gradients associated with field fluctuations cause the density 
perturbation to oscillate and decay. Also shown are the dark en- 
ergy equation of state wq = Pq/pq and the adiabatic sound speed 
squared c^ = Pq/pq (not to be confused with the effective sound 
speed c e = 1) which indicate the field is Hubble drag dominated 
until recently (see Appendix). 



and hence from Eqn. JSJl 
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where the two pieces are the adiabatic term and the 
isocurvature term. The integrals in Eqn. (|20|l can be 
expressed in terms of hypergeometric functions. In the 
fiducial model 
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Moreover since the field only experiences a small range 
in the potential throughout the whole expansion history, 
any underlying form of Vq that satisfies these require- 
ments will have the same phenomenology. We find that 
in the context of the fiducial model €q < 0.6 is required 
for«) Q (o=l)<-2/3. 

Given potential energy domination, the energy and 
pressure fluctuations are related to the field fluctuations 
as 



= 5qPq = -$Pq = Vqi 
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and remain nearly constant during the evolution. The 
superhorizon evolution of the comoving curvature from 
Eqn. (J7J is given by 
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Combined with the requirement that A<I> s=s — Q this re- 
lation implies that we should set the initial dark energy 
fluctuations to be 5qi ~ —8Q to cancel the quadrupole 
in the fiducial model. 

While this condition is roughly correct, the scales 
that are responsible for the quadrupole in the SW ef- 
fect (k ~ 0.0002 Mpc -1 ) are on the horizon scale today. 
Since the quintessence field has a sound horizon equal to 
the horizon, the fluctuations in these modes will have al- 
ready begun to decay from their initial values. In Fig. 0] 
we show the time evolution of the dark energy density 
perturbation. Because the argument above might appear 
to require wq = — 1 exactly, we have chosen to illustrate 
the behavior in a model with wq(o, = 1) = —0.94 and 
£q = 0.18 with the other parameters equal to their fidu- 
cial values. Since wq — ► — 1 rapidly with redshift, the 
change in £>* from the fiducial model is only 0.4%. As 
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FIG. 5: Quadrupole transfer functions in the fiducial adiabatic 
model ("A") and a model with additional isocurvature perturba- 
tions of Si = -15 ("A&I"). The isocurvature ISW effect (ilSW) 
cancel the SW effect for the temperature quadrupole from large 
scales while leaving the polarization nearly unchanged. 



the background evolution is nearly indistinguishable from 
the true fiducial model, we will employ this choice for the 
isocurvature analog of the fiducial model. 



In practice we have taken a potential Vq 



'Q l 



72 



with mq = 10~ 42 GeV and an initial position consistent 
with €q and dQ/d\na = initially. For scales near the 
peak of the SW effect in the quadrupole the density per- 
turbation has decayed by about a factor of two by z ~ 1. 
Consequently an initial density perturbation of 
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should be optimal for reducing the SW quadrupole. We 
call models with this type of fully correlated adiabatic 
and isocurvature models "A&I" models. 

In Fig. |31 we show the effect of this initial condition 
on the gravitational potential. In this case the change in 
the gravitational potential A<& rs — Q and one achieves 
the desired effect of eliminating the SW effect in the 
quadrupole. 

The reduction of the SW quadrupole must not come at 
the expense of an enhancement in the ISW contributions 
to the quadrupole. Fortunately, this is a natural con- 
sequence of the effective sound speed of the scalar field 
c e = 1. Scales near the peak contribution of the adiabatic 
ISW effect are well within the horizon by dark energy 
domination. Consequently as can be seen in Fig. 0] any 
initial isocurvature perturbation in the dark energy will 
have decayed before dark energy domination. In Fig. |3 
we show the evolution of the potential given the initial 
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FIG. 6: Temperature power spectrum of the A&I model with sev- 
eral choices for the initial isocurvature amplitude compared with 
the WMAP data. The shaded region corresponds to the 68% and 
95% cosmic variance confidence region of the fiducial A&I model 
with Si = —15 (solid line). 



conditions of Eqn. IJ22JI . Note that despite the large ini- 
tial isocurvature perturbation it has almost no effect on 
the potential evolution and hence the contributions to 
the quadrupole. 

In Fig. [SI we show the temperature and polarization 
transfer functions of the isocurvature model. The isocur- 
vature ISW effect almost perfectly cancels the SW effect 
on large scales leaving a quadrupole that consists almost 
solely of the adiabatic ISW effect. Note that the isocurva- 
ture conditions leave the polarization essentially unmod- 
ified as expected since they only change the potentials at 
z < 2 in Fig. 

In Fig. we show the temperature power spectrum of 
this model for several choices of Si with all other param- 
eters held fixed. Note that for Si ~ —15 the suppression 
is fairly sharp around the quadrupole in spite of scale 
invariance in both the adiabatic and isocurvature initial 
conditions. The polarization auto (EE) and cross (QE) 
spectra are shown in Fig. [7| Notice that although the 
EE spectrum remains nearly unchanged, the QE spec- 
trum has a reduced cross correlation with a sign change 
at the quadrupole. (Note that our sign convention for E 
is opposite to CMBFAST.) Because the cross spectrum 
in the adiabatic case arises from the correlation between 
the SW quadrupole and the polarization, it is affected by 
the cancellation of the SW quadrupole as the product of 
the two transfer functions in Fig. [5] show. 

It is interesting to compare the isocurvature method 
of lowering the quadrupole to other possible solutions. A 
related mechanism involves modifying the sound speed 
of the dark energy [37], LlU HE 03 ■ Here one retains adi- 
abatic initial conditions but modifies the effective sound 
speed of the dark energy (see Appendix) . The dark en- 
ergy then contributes density fluctuations between the 
horizon and the sound horizon and reduces the adiabatic 
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FIG. 7: The polarization EE and cross correlation QE spectra 
for the fiducial adiabatic model compared with the fiducial A&I 
model of Fig. 151 with Si = —15. Although the EE spectrum re- 
mains essentially unchanged, the correlation is modified at the low 
multipoles due to the reduction of the large scale contributions to 
the quadrupole from the SW effect which are the source of the 
correlation in the adiabatic model. 



FIG. 8: Alternate models for the low quadrupole. The adiabatic 
ISW effect may be reduced by lowering the effective sound speed 
of the dark energy, here to c e = 0.03. The SW effect can be largely 
eliminated through a cut off factor C(k) in the initial adiabatic 
power spectrum, here chosen to remove power for k < fe cu t = 0.0005 
Mpc . Although the latter does not modify the transfer function, 
we have illustrated its effects by showing T(C 1 ' 2 (k). Note that the 
cut off affects the polarization whereas the sound speed does not. 



ISW effect by suppressing the decay of the gravitational 
potential. In Fig. [S] we show the transfer functions for 
a model with wq = —2/3, c e = 0.03 with Qq — 0.37 
and tt m h 2 and n^h 2 held fixed to the fiducial values such 
that -D* and the shape of the peaks remain the same as 
in the fiducial model. Note that unlike in the isocurva- 
ture case, the SW contributions to the quadrupole are 
unmodified but the adiabatic ISW contributions are re- 
duced. The resulting effect on the temperature power 
spectrum is shown in Fig. |5J However, like the isocurva- 
ture case, the polarization is largely unchanged since the 
dark energy again only affects low redshifts. A qualita- 
tive difference appears in the cross power spectra, which 
remains unchanged in this case (see Fig. llO|) . 

Finally, the quadrupole can be lowered by removing 
power on scales associated with the SW effect. Here we 
take the model I23| 



A^(fc)=A^(fc)| fid C(fc), 

C(k) =1 - e -(*A°ut)"-' 



(23) 



with fc cut = 0.0005 Mpc^ 1 and n cut = 3.35. The trans- 
fer functions for this model are the same as the fidu- 
cial model but for illustrative effect we plot them as 
Ti(k)C l / 2 (k) in Fig. Like the isocurvature model, the 
quadrupole in Figgis suppressed due to the elimination 
of the SW quadrupole. Unlike the isocurvature case, both 
the EE and QE spectra are suppressed since the tem- 
perature quadrupoles are also absent during reionization 
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FIG. 9: Temperature power spectra of the sound speed and cut 
off models of Fig. The sound speed model lowers the adiabatic 
ISW contributions to the quadrupole. The cut off model reduces 
the SW contributions to the quadrupole. 



(seeFig.HTJ. 

IV. LIKELIHOOD ANALYSIS 

In this section we check whether the adiabatic plus 
isocurvature (A&I) model of the previous section is fa- 




FIG. 10: Polarization spectra for the sound speed and cut off mod- 
els of Fig. El The sound speed model changes both spectra negli- 
gibly whereas the cut off model suppresses power in both spectra 
simultaneously. 



TABLE I: Scale Invariant ACDM Model Parameters: 
60 & 0E Data. 



WMAP 



Parameter 



Maximum Likelihood 



Matter / Critical Density, Q m 0.27 

Baryon / Critical Density, fif, 0.046 

Hubble Constant, h 0.72 

Optical Depth, r 0.17 

Initial Amplitude, 5 (i 5.07x 10" 5 



TABLE II: Best fit (assuming scale invariance) \ 2 an< l degrees of 
freedom (DOF) values using WMAP data. 



Model 


Data 


x 2 


DOF 


Fiducial (Table [jj 


ee 


976 


894 


c e = 0.03 


ee 


973 


893 


Cut-off (fc cut = 0.0005Mpc" 


- 1 ) ee 


972 


893 


Adiabatic & Isocurvature 


ee 


972 


893 



Fiducial (Table 

c e = 0.03 

Cut-off (fc cu t = O.OOOSMpc- 1 ; 

Adiabatic & Isocurvature 



66 + QE 1430 1343 

ee + QE 1428 1342 

ee + QE 1428 1342 

ee + QE 1427 1342 
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FIG. 11: Areas enclosing 95% of the probability for the dark 
energy isocurvature perturbation (Si) and the optical depth (t). 
The black area is using only the ©0 WMAP data and the gray 
area uses both the 00 and BE WMAP data. 



vored by the WMAP data [l|. For simplicity we will 
assume a scale invariant spectrum, zero background cur- 
vature and zero neutrino mass. Then the maximum like- 
lihood fit with no dark energy perturbations is given in 
Table [TJ When the spectral index is fixed the only pa- 
rameter with enough freedom to significantly effect the 
low multipoles is the optical depth, r. 

Keeping the other parameters in Tabled fixed we vary 
the magnitude Si of the primordial isocurvature pertur- 
bation in dark energy relative to the comoving curva- 
ture between —60 < Si < 25 and the optical depth 
< t < 0.29 on a 28 x 30 grid. The amplitude in Table|I| 
is scaled by exp(r — 0.17); recall that the power spectrum 
is scaled as the square of this quantity [see Eqn. ©]. 
This is to take into account the well known degeneracy 
between r and the amplitude. For each grid point the 
likelihood was evaluated using the software provided by 
WMAP [42, 13 and the resulting two dimensional sur- 
face was interpolated. A uniform prior in both parame- 
ters was assumed. 



The maximum likelihood points are Si = —15.3 and 
t = 0.16 for WMAP 69 data only and Si = -12.1 and 
t = 0.18 for WMAP 96 and QE data. The x 2 for 
these parameter values and other models are giving in 
Table ITU The reduction in \ 2 is achieved by only modi- 
fying the first few multipoles of the predicted spectrum. 
Areas enclosing 95% of the probability with and without 
the QE data are shown in Fig. ^] The area perime- 
ters lie on constant probability contours. As can be seen 
from Fig. ^2 the weight of the probability distribution 
favors Si < 0, i.e. an isocurvature density perturbation 
that is negatively correlated with the comoving curva- 
ture perturbation, Q. This preference comes from the 
fact that the lowered temperature quadrupole in such a 
model better fits the observations. The addition of the 
BE data reduces the favored magnitude of the isocurva- 
ture perturbation Si and slightly increases that of r as 
QE is suppressed in the model but not in the data |49J |. 
Note that the restriction to scale invariant models does 
not affect the relative improvement that Si makes for a 







/'\00+0E 


0.05 






0.04 


00// 




0.03 






0.02 




\\ 


0.01 


/ / 






i i , ,i 



-60 



-40 



-20 



20 



FIG. 12: Probability density functions for the dark energy isocur- 
vature perturbation using WMAP 00 (solid) and WMAP 00 and 
QE data (dashed). 



given r since across the scales relevant to the first few 
multipoles any nearly scale invariant model will have the 
same shape. Allowing the tilt to vary does allow low r 
models if only the 00 data is employed but these are 
eliminated with the inclusion of QE data (c.f. |2j and 

One dimensional probability distributions obtained by 
integrating out r are displayed in Fig. 1121 The mean and 
68% confidence regions for Si are -15.2 ± 8.2 for WMAP 
00 data only and -11.8 ± 7.1 for WMAP 00 and QE 
data. The probability of Si > is 0.04 using WMAP 
00 data only and 0.06 using WMAP 00 and QE data. 
This is in sharp contrast to the estimates of the amount 
of totally correlated CDM, baryon or neutrino isocur- 
vature perturbation as estimated in |5li l52j where the 
probability for a positive isocurvature perturbation was 
roughly 50%, i.e. the isocurvature perturbation was not 
significantly different from zero. The reason for this is 
that correlated non-dark energy isocurvature perturba- 
tions have a similar effect on a much larger range of low 
£ multipoles as they rely on a simple reduction of the 
SW effect. The dark energy introduces a scale associated 
with the current horizon and so the dark energy isocur- 
vature ISW effect can only reduce the SW effect on the 
large scales associated with the quadrupole. 



V. INFLATIONARY CONTEXT 

As we have shown, it is essential that the dark energy 
isocurvature perturbation is anti-correlated with the adi- 
abatic perturbation so that its ISW effect can coherently 
cancel the SW effect. Correlated or anti-correlated adia- 
batic and isocurvature perturbations result when the adi- 
abatic curvature perturbation itself is generated by the 
isocurvature perturbation J53, |54j. Ref. J3g suggested 
a three field model based on the curvaton mechanism 



|55l |56j, |57j in order to produce anti-correlated isocurva- 
ture dark energy perturbations. 

Another way of establishing the anticorrelation is 
through the variable decay mechanism |58l [59j where 
there is a second light scalar field during inflation which 
determines the decay rate T of the inflaton. In this 
method, the resulting comoving curvature perturbation 
is given by 






1ST 

6 r 



(24) 



If we associate the variable decay field with the 
quintessence field, Q, and take T oa Q then = 
— (1/6)SQ/Q. Then, if we further take the quintessence 
potential to be Vq oc Q 2 , we have §Q i = 25Q/Q. There- 
fore 



Sq, = -120 



(25) 



which was the mean value found in iJIVI when both 
WMAP 00 and BE data were used. 

Unfortunately, a problem with this and other inflation- 
ary scenarios arises because of the gravitational waves 
generated during inflation. Both the scalar field and 
gravitational waves are nearly massless degrees of free- 
dom during inflation and hence acquire closely related 
quantum fluctuations. The contribution of th e g ravita- 
tional waves to the 00 spectrum for £ ^ 1 is [6(| 



l(l + l) r ee 



2irG 



H-, 



inf 



2vr 



(26) 



where _ffi n f is the Hubble parameter during inflation. On 
the other hand, Eqn. Q implies that the curvature fluc- 
tuations induce a SW effect of 



^+l)^ee 



2tt 



-a 



SW 



25 Cl 



(27) 



Now A? is related to the inflationary power spectrum of 



A 2 - 
A <5Q - 



H; 



inf 



2tt 



(28) 



through 



120 



V 6 



= ^WttGcqSQ , 



(29) 



where cq was defined in Eqn. I|15|) . Unlike inflaton cur- 
vature fluctuations, these curvature fluctuations are sup- 
pressed with a small cq because the cancellation condi- 
tion involves the energy density fluctuation in Q. There- 
fore 



1(1+1) 

2tt 



a 



s\v 



G-k ( H inI 



225 



2tt 



(30) 
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The predicted ratio of gravitational wave to SW contri- 
butions is r ft 450/eQ. 

As discussed in illlll an equation of state today of 
wq ft —1 requires a small slow roll parameter for the 
quintessence (cq < 0.6) at the field position today. The 
current observational limit is about r < 1 |2|< There- 
fore, the isocurvature perturbations in a nearly frozen 
field model for the quintessence field could not have orig- 
inated from this inflationary mechanism. 

Note that the problem of excess gravitational waves is 
independent of the precise potential or inflaton decay rate 
relation as it arises from the phenomenologically required 
condition Eqn. IJ25)) . It applies to any model (including 
the curvaton model proposed in Ref. J36J) for which the 
field fluctuations SQ at dark energy domination are di- 
rectly associated with the inflationary fluctuations of a 
light canonical scalar field [Eqn. (|28[) ]. 

The only way to evade this conclusion is for SQ to grow 
by a factor of 30 or more between horizon crossing during 
inflation and today. It is not clear how to do this for a 
canonical scalar field whose perturbations were generated 
during inflation. 



VI. DISCUSSION 

We have explicated the mechanism by which scale in- 
variant isocurvature perturbations in the dark energy can 
lead to a sharp suppression of the quadrupole tempera- 
ture anisotropy [3a| ■ The Integrated Sachs- Wolfe (ISW) 
effect in the dark energy dominated epoch can coherently 
cancel the Sachs- Wolfe (SW) effect from recombination if 
its energy density fluctuations are set to be strongly an- 
ticorrelated with the initial curvature fluctuations. We 
call these A&I models. Assuming a scale invariant spec- 
trum and flat background, we detect the presence of A&I 
perturbations at the 95% confidence level. 

As is well known [6l|, |62, |63 (see also Appendix), 
isocurvature modes in the dark energy rapidly decay for 
tracking models. Models in which the scalar field and 
its isocurvature perturbations are essentially frozen |62( 
due to a shallow potential slope relative to the Hubble 
parameter are better suited. We have shown that the 
requirements on such a model comes mainly through the 
quintessence "slow roll" parameter cq. 

As for the origin of the A&I correlation, we can asso- 
ciate the dark energy with a variable decay rate of the 
inflaton. We find that the right level of isocurvature per- 
turbations is naturally predicted. Unfortunately, the re- 
quirement of a shallow slope of the potential or small cq 
implies too high a level of gravitational waves. For this 
model to work some mechanism is required for amplifying 
the field fluctuations by a factor of 30 or more between 
inflation and dark energy domination. 

A&I models should be contrasted with those that in- 
troduce a cut off scale to the perturbations. A conceptual 
problem of the latter class is that it introduces the sharp 
reduction at the quadrupole "by hand" . In other words, 



there is a new coincidence problem between the cut-off 
or topology scale and the horizon today. There is also no 
significant evidence for features at any other scale (e.g. 
|64J |65j). Phenomenologically, by eliminating the per- 
turbations altogether, the cut off models also eliminate 
the source of large angle polarization. Thus, the cut-off 
and A&I models are potentially distinguishable from the 
polarization autocorrelation (or EE) spectrum. On the 
other hand both models predict a reduced temperature 
polarization cross correlation. 

In principle, these alternatives can also be distin- 
guished by large scale measurements of the density field 
as a function of redshift and its ISW correlation with the 
CMBjplElMIe^, e.g. through highredshift galaxies 
HI ETIzlLZill or cosmic shear [Zilzl |?|, but mea- 
suring the small signals involved will require exquisite 
control over systematics in the surveys. 

Both cut off and A&I models operate by reducing 
the SW contributions to the quadrupole. Unfortunately, 
the SW effect only contributes approximately half of the 
quadrupole in adiabatic models with a cosmological con- 
stant. The remaining portion comes from the ISW effect 
and receives contributions across a wide range of subhori- 
zon scales. Neither A&I nor cut off models can suppress 
the ISW effect at the quadrupole. In the former, the 
high sound speed of the dark energy prevents substantial 
density perturbations on subhorizon scales. In the latter, 
a cut-off on small enough scales to affect the ISW effect 
would remove too much small scale power and distort the 
higher £ multipoles. 

The adiabatic ISW effect can be modified by changing 
the dynamics of the dark energy by lowering its effective 
sound speed. Alone it only amounts in a fairly small 
reduction if the other multipoles are not to be adversely 
affected. This model is also distinguished from the A&I 
and cut off models in that it affects neither the auto nor 
the cross correlation spectra of the polarization. This 
is because it is the SW quadrupole that is responsible 
for the polarization and hence its correlation with the 
temperature. 

In combination with an anticorrelated A&I model, a 
sound speed modification alters the multipole at which 
the SW effect is canceled. We find that the canonical 
sound speed c e = 1 is nearly optimal in producing a sharp 
suppression at the quadrupole though a slight increase in 
the sound speed can actually lead to a slightly sharper 
reduction. 

Since the first COBE detection, the low quadrupole 
temperature anisotropy in the CMB has provided a tan- 
talizing hint that new physics may be hovering on the 
horizon scale. With upcoming polarization auto and 
cross correlation data from WMAP, we may soon more 
than double the information on this intriguing problem. 
All of the alternatives discussed here have distinct, albeit 
cosmic variance limited, predictions for these spectra. 
The predictions are especially distinct at the quadrupole 
and octopole but it remains to be seen how well these 
large-angle polarization fluctuations can be separated 
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from galactic foregrounds and instrumental effects. 
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APPENDIX A: DARK ENERGY 
PERTURBATIONS 

1. Covariant Conservation 



2. Equation of State 



Following |43j, [77J , let us define the most general per- 
turbation to the Friedmann Robertson Walker (FRW) 
metric as 



g°" = -a- 2 (l-2A), 
g - -a B , 



9 



a-\^ -2H Ll l - 



21*30, 



(Al) 



rpO 

1 


= ~PQ, 




i 


= QQi, 




3 


= PQ &t 3 


+ n Q i j 



where 7^ is the FRW 3-metric of constant (comoving) 
curvature K = Hq(Qt — 1)- Likewise let us define the 
most general stress-energy parameterization of an energy 
density component Q 



(A2) 



where pq is the energy density, qQi is the momentum 
density, pq is the pressure or isotropic stress and IIq 1 is 
the anisotropic stress of Q. The anisotropic stress is de- 
fined as the trace free portion of the stress such that 
IIq 1 ~0. If Q is to dominate the expansion at any time 
then qQ and Hq must vanish in the background due to 
isotropy. 

If Q is non-interacting or "dark" then this stress energy 
tensor is covariantly conserved. Conservation then leads 
to general equations of motion for the stress-energy com- 
ponents. Retaining terms linear in the metric fluctua- 
tions where they combine with stress energy components 
(pq,Pq) in the background, we can reduce VT^, = to 



d h 

dr\ a 



' d 


.d 




+ 4- 


di] 


a 



PQ 



QQi 



3-pq - V^ Q 
(p Q +P Q KV t B l + 3H L ), 

(PQ+PQ)ViA. (A3) 



The conservation equations represent a general but in- 
complete description of the dark energy as they leave the 
spatial stresses pq and IIq unspecified. 



The isotropic stresses may be rewritten in terms of an 
equation of state 



Pq(v, Xi) = wq(v> %i)pQ(v> %i) , 



(A4) 



where wq is in general a function of position and time. 

For the scalar degrees of freedom it is useful to recast 
Eqn. (|A3|) in terms of an expansion of the fluctuations 
into scalar harmonic modes defined by the complete set 
of eigenfuctions Y of the Laplace operator 



V 2 F 



-k 2 Y. 



(A5) 



The spatial fluctuations in each mode are given by four 
mode amplitudes Sq,5wq,uq,itq as 

$PQ =SqYp q , 

Sp Q =(wq5q + 5w Q )Yp Q , 



QQi =u Q (~k 1 V i Y)p Q , 



n 



Q 3 



E7T Q (£f 



'v'Vj-y - 



l^y)pQ 



(A6) 



Note that for a spatially flat metric Y = e lkx and the 
mode amplitudes are the Fourier coefficients of the fields. 
In the literature (e.g. [43J]'). one often defines instead a 
quantity Vq, related to the bulk velocity or energy flux, 
such that uq = (1 + wq)(Vq — B). In Eqn. l|A6p and 
throughout the remainder of this section /?q, pq, and 
wq = Pq/pq are to be reinterpreted as the average or 
background quantities and depend only on time. Here- 
after we will typically omit the spatial harmonics Y by 
assuming a harmonic space representation of the pertur- 
bations. 

The conservation equations then become 

Sq=— 3-5w q - ku Q - (1 + wq)(kB + 3H L ) , 

uq =-(3wq - l)u Q + k(w Q 5Q + Swq) 

-^w Q (l-3K/k 2 )kn Q + (l + w Q )kA, (A7) 

where we have used the conservation of energy relation 
in the background [see Eqn. (|A3|) ] 



PQ = -3-(l + Wq) P q. 



(A8) 



Here we have followed the same convention for the 
harmonic representation of the metric fluctuations, e.g. 
A(r},Xi) = A(T],k)Y. 

The quantity Swq represents spatial fluctuations in the 
equation of state and specifies the pressure fluctuation as 



Sp Q = (w q 6q + 5wq)pq 



(A9) 



Equation of state fluctuations can arise from intrinsic 
internal degrees of freedom in the dark energy or from 
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temporal variations in the background equation of state 
through a general coordinate or gauge transformation. 

The gauge transformation is defined as a perturbation 
in the temporal and spatial coordinates of x^ = a; M + 
(T, LY l ) with respect to the background under which 





.4: 

B 


= A-T - -T, 
a 

= B + L + kT, 




H L 


= H L 


- k -L- h -T, 
3 a 




Ht - 


= ix J 1 


+ kL, 


the metric 


and 








SpQ = 


-' S Pq 


- PQ T , 




5pq = 


-' S PQ 


-PqT, 



(A10) 



U Q = UQ - k(l + w Q )T , 



(All) 



for the dark energy variables. 

Starting from the uniform density gauge where there is 
no dark energy density fluctuation SpQ — 0, there arises a 
contribution to the pressure fluctuation under the gauge 
transformation [Eqn. (|A11|) ] of 



$Pq - S PQ = SpQcl ■ 
where the adiabatic sound speed is 

2 = PQ_ _ dp Q /d\na 
a p Q dpQ/dlna 



(A12) 



=W Q - o 



ldln(l 
3 



WQ) 



d\na 



(A13) 



If in the uniform density gauge the pressure fluctuation 
also vanishes (SpQ = 0), the dark energy only has adi- 
abatic pressure fluctuations and from Eqns. (|A9() and 
(IAT21) 



SWQ adi =(c 2 a ~Wq) 



(A14) 



In this case, the pressure fluctuation and gravitational 
terms in the momentum conservation equation (|A7|) be- 
come 



kc 2 a S Q + (l + w Q )kA. 



(A15) 



If the dark energy accelerates the expansion wq < —1/3. 
If wq is also slowly varying, the adiabatic sound speed 
is imaginary: c 2 ~ wq < 0. In the small scale New- 
tonian approximation (and gauge), the Poisson equation 
gives A ~ —SQ/(krj) 2 during the dark energy dominated 
epoch. Hence for c 2 < pressure gradients would en- 
hance potential gradients in generating momentum den- 
sity and the dark energy would collapse faster than the 
dark matter. Thus for gravitational instability in the 
dark energy to be stabilized by pressure gradients a 
source of non-adiabatic pressure fluctuations is required. 



These can be supplied by internal degrees of freedom in 
the dark energy. 

Following |ig , the role of the background equation of 
state parameter as a closure relation between the pres- 
sure and density can be generalized for an inhomogeneous 
dark energy component. If one requires that the pressure 
fluctuation is linear in the energy and momentum density 
fluctuations, general covariance [or gauge invariance for 
linear fluctuations, see Eqn. (|All|l ] requires that 



<5wQ|non-adi <* Sq +3-— . 

a k 



(A16) 



We choose to specify the proportionality through an ef- 
fective sound speed such that the full equation of state 
fluctuation is given as 



Sw Q =( c a - w Q )S Q + (c\ - c 2 a ) [Sq + 3 



■UQ 



(A17) 



Note that c 2 , — c 2 for adiabatic pressure fluctuations. 
The quantity c 2 may in principle be a function of time 
and k. 

With this ansatz for Swq, the covariant conservation 
equations (|A7ll become 



Sq =-3-(c e -wq)Sq 



~) (<*-«&¥ 



ku Q -(l + w Q )(kB + 3H L ), 



: [3(wq 



2 2\ 

c: - c„) 



1]uq + kc 2 e Sc, 



- 3 w Q (l - 3K/k 2 )kir Q + (1 + w Q )kA . (A18) 

Note that unlike in Eqn. 1|A15(1 . it is c 2 and not c 2 that 
is the coefficient for Sq in the momentum conservation 
equation; hence c 2 controls the stability of density fluc- 
tuations. More formally in the zero momentum gauge or 
rest frame of the dark energy where uq = 0, Eqns. (|A9|) 
and l|A17|l imply 



Swr 



(c e - WQ)S Qrcst , 

rest — 



c 2 e Sp Qrcst 



(A19) 



so that c 2 is the sound speed in the rest frame. This 
property is useful in the calculation of c 2 for a given dark 
energy model. 

Likewise, a closure relation between the anisotropic 
stress 7tq and the energy and momentum density can also 
be supplied through a generalized equation of state [46| . 
However unlike the non-adiabatic stress, a non- vanishing 
value is not required for stability by wq < —1/3. We 
shall see now that for scalar field dark energy candidates 
it in fact does vanish in linear theory. 
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3. Quintessence and k-Essence Fields 

The stress-energy tensor for a scalar field Q with a 
Lagrangian F(X, Q) is J23 



Tt = |^VQV v Q + F*t ; 



where the kinetic term 
X = 



1. 



rV^Qv.g. 



(A20) 



(A21) 



For a scalar field with the canonical kinetic term in 
Eqn. (|A22() . cjf = 1. More generally, to achieve a con- 
stant c 2 one requires 



dF 
dX 



X \~^ 



(A29) 



where Vr- is a constant with dimensions of energy density. 
If for simplicity one further assumes that the kinetic and 
potential terms are additive dF/dQ = —OVq/SQ and 



For the case of a canonical kinetic energy term 
F(X,Q) = X-V Q , 



(A22) 



where Vq(Q) is the scalar field potential. Q is then called 
the quintessence field. For the case F(X, Q) = —X — Vq, 
the dark energy is called a phantom field (e.g. |jj|) an( i 
for a more general modification of the kinetic term, a 
k-essence field |80) . 

Matching terms between Eqn. (|A2(I and l|A20j) . and 
dropping contributions that are second order in the spa- 
tial gradients of the field 

dF 
PQ = F, 



IV, = o . 



OF 



(A23) 



In particular the density and pressure fluctuations be- 
come 



, d 2 F dF, r 

dF dF 

5p *=8X 6X+ 0Q dQ > 



d 2 F dF\ . 

2 dQdX X ~dQ) Q ' 

(A24) 



where 



6X = a- 2 (Q5Q- AQ 2 



(A25) 



Note that under a gauge transformation Q transforms 
as a scalar 



and 



5Q = 6Q-QT, 



dF _ 2 • 
Pquq = -^rza QkSQ . 



(A26) 



(A27) 



Hence the uniform field gauge where SQ — coincides 
with the rest gauge uq = 0. Since this is the gauge in 
which c 2 = Spq/Spq, the effective sound speed of a scalar 
field is El, 



dF/dX 



2(d 2 F/dX 2 )X + (dF/dX) 



(A28) 



F = 



2< 2 



1 



v k -v Q 



(A30) 



The effective sound speed of Eqn. (|A28|I along with 
wq in the background closes the evolution equations 
Eqn. l|A18Jl once the metric fluctuations are specified by 
the Einstein equations (see AA50 . The evolution of dark 
energy density perturbations is then completely defined 
by a choice of initial conditions for Sq and uq. 



4. Background Evolution 

To obtain the background evolution of the dark energy, 
one begins with the energy density conservation equation 
in terms of the field variables 



Q + (3c;; 



n-.. + ^(n 



a 

)-( 
a 



_1 dp Q 
dQ 



or 



d 2 Q 
d In a 2 



dQ 



a^-^^te + ^l- 



£1 

H 2 



dF 
dX 



(A31) 



= 0. 



For the canonical kinetic term, 

tPQ 3 dQ 

dlna 2 2 l WT) dlna~ 



1 dV Q 
H 2 dQ 



0. 



(A32) 



(A33) 



where wt = Pt/ ' Pt is the equation of state for the sum of 
all components. Note that wt depends on Q and hence 
in practice it is more convenient to solve a set of coupled 
first order differential equations in Q and HdQ/dlaa = 
y2X. Likewise, since 1 + wq appears in the evolution 
equations, one calculates this quantity directly as 



1+WQ 



2{dF/dX)X 
2{dF/dX)X - F 



(A34) 



Combined with Eqn. I|A29(1 . this relation implies that if 
the field is potential energy dominated wq — > — 1 and if 
it is kinetic energy dominated wq —> c 2 . 

With equation (|A32|I . the expression for the adiabatic 
sound speed becomes 



2 c 2 e dpQ/dQ - dF/dQ 

C p -f" 



3H 2 {dF/dX){dQ/dlna) 



(A35) 
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The adiabatic sound speed takes on a specific form in 
the case where the Hubble drag inhibits the motion of 
the field. In this case Opq/OQ is nearly constant. If wt 
is constant, dQ/d\na reaches a "terminal velocity" that 
is independent of the initial conditions and scales with 
H . Utilizing the constant effective sound speed form for 
F in equation i|A30(l , the equation of motion l|A32|) can 
be written in the form of the Bernoulli equation and the 
adiabatic sound speed becomes 



cl = cl--(c 2 e + l)(3 + w T ). 



(A36) 



For the case of a canonical kinetic term c 2 a = — 2 — wt 
H2 (see also Fig. EJ. 

On the other hand, if the Hubble drag is negligible 



dQ 
din a 



-3(2cl-l-w T )/2 



(A37) 



and assuming kinetic energy domination wq — c 2 , = 



5. Initial Conditions and Evolution 

We numerically solve the dark energy evolution equa- 
tions in comoving gauge where the total momentum den- 
sity vanishes 



Ptut 



, PiU., 



= 0. 



(A38) 



Here i runs over all species of energy density. We define 
the total density and pressure fluctuation similarly. The 
auxiliary condition Ht = completely fixes the coordi- 
nate freedom. We call the remaining metric degrees of 
freedom ( = Hl, £ = A and Vt = B, where Vt also has 
the interpretation of the total momentum weighted ve- 
locity [see J43J and the discussion below Eqn. I|A6() ] . The 
Einstein equations become 



c- 



K 


a 






a 


1 




a pr +Pt 


V t 


=(*- 





-Sp T 



3K 

T2- I PT^T 



(A39) 



where recall K is the background spatial curvature. The 
continuity equation for the total density perturbation 
completes the basic equations 

S T = -3-Sw T - (1 + w T )(kV T + 30 , (A40) 

a 

as the momentum conservation equation was used to re- 
late £ to the stresses in Eqn. (|A39(I . In that equation 



(A41) 



(fc 2 - 3iT)$ = 4irGa 2 p T 6 T 



and is equal to the curvature fluctuation in the Newto- 
nian or longitudinal gauge $ = -Hz, | Newt- For reference, 
the Newtonian potential \& = A|Newt employed in JiTlis 

* = -$ - 87rGa 2 jj T 7r T /fc 2 , (A42) 

and the total velocity obeys an auxiliary relation 



V T + -V T = -k(£ - *) 
a 



(A43) 



that is useful for checking numerical solutions. 

Equation l|A39|) implies that the curvature fluctuation 
C changes only in response to stress fluctuations. As in 
the discussion of dark energy perturbations above, the 
evolution equations are closed through an assumption 
for the stress perturbations. In practice, for numerical 
solutions that extend to subhorizon scales, we replace 
Eqn. (|A40p with the conservation equations for the indi- 
vidual energy density species to implicitly solve for 5wt- 

For the initial conditions, assuming nearly adiabatic 
stresses 



fipr ; — Spr 

Pt 



< \6pi 



(A44) 



and a closure relation for the anisotropic stress of the 
neutrinos which comes from the Boltzmann equation [e.g. 
[H Eqn. (12)] 



n v = jkriVr* 



(A45) 



we obtain the solution to the evolution equations in the 
radiation dominated regime 



4 1 + 3oy K " 2 

* T= 9iT2^ (1 - 3 F )(b?)Ci ' 



V T =-l: ' 



3 1 + 2a„ 



(hri)(i , 



(A46) 



where a„ = 2p v jVz>p r accounts for the anisotropic stress 
of the neutrinos 



Pv 

7TT = 7IV 

Pr 



2a, 



l + 2a, 



-(M 2 G 



(A47) 



Here p r = p^ + p v , the total radiation density and 
we have kept only leading order terms in power of krj. 
Given that the stress perturbation Spt/(pt + Pt) — 
0[(kri) 2 ]Ci, Eqn. (|A39|I shows that the curvature per- 
turbation is nearly constant for superhorizon adiabatic 
stresses C{k,Vi) = Q(k) [43j - 

The dark energy perturbations associated with the cur- 
vature fluctuation can then be related to the total density 
perturbation by substituting these relations back into the 
Einstein equations i|A39D for £ and £ and employing them 
in the conservation equations l|A18(l . The result is 



Sq = A s (l + wq)5 t , 
u Q = A u (1 + WQ){kri)6 T 



(A48) 
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where assuming a radiation dominated initial condition 



A u — 



3 C 2 + (l + 6a l/ )(3w Q -2) 



A s 



4(1 + 3a„)[8 + 3c2(2 - 3c 2 a + 3wq) 

3 1 - 6(cg - c 2 a )A u 

4 1 - 3(w Q - C 2)/2 • 



12w Q ] ' 
(A49) 



We have here assumed that wq and c\ are nearly constant 
compared with the expansion time but have allowed (1 + 
wq) to vary as appropriate for a scalar field under the 
Hubble drag by factoring this quantity out of Eqn. I|A48|) . 

The curvature initial conditions for the dark energy 
take on this rather intricate form involving A u in the 
relation for Ag since the internal non-adiabatic stress of 
the dark energy cannot be neglected even though its con- 
tribution to the total non-adiabatic stress in Eqn. (|A44|) 
can be ignored. In general, the curvature mode always 
carries non-adiabatic stresses beyond the leading order 
in (krj) and pi/p r where i ^ 7, v. Formally these vanish 
if the initial conditions are taken to r\i — » 0. 

Deviations in the initial conditions for the dark en- 
ergy from these relations represents an isocurvature mode 
since the dark energy is assumed to carry a negligible 
fraction of the net energy density at the initial condi- 
tions so that the radiation density fluctuation S r = St- 
If we also assume that the relative deviations are large 
5Q(r]i,k) ^ (1 + WQ)ST(i]ijk) then the dark energy- 
radiation entropy fluctuation 



(1 + wq)Sqt =S q - — — Q-8 r , 

1 + W r 



(A50) 



which corresponds to the usual definition of the isocur- 
vature mode as being generated by SQ r . The distinction 
here is that the adiabatic mode has Sqt — 0(6t) 7^ 
due to the intrinsic entropy of Q and evolution of wq . 

Neglecting the metric fluctuations generated by the 
dark energy fluctuations, we find that the the evolution 
equations (|A18|I are solved by a linear combination of the 
adiabatic mode and 



where 





6 Q = J 8 (knY, 






u Q = i u (k V r +i , 




h 


[■ ! 3(1 + w T + 2c 2 a - 2c\ - 


-^wq)! 


In 


1 + 3wt 





and p solves the equation 



G 



1 + 3wt 



(Ce " W Q ) 



In 



1 + 3u>t 



(A51) 



(A52) 



(c2 - cl) . 
(A53) 



For initial conditions in the radiation dominated era, the 
solutions to this quadratic equation are 



P = -1 - 2 C » + Sw Q ± VO- + 3c a/2) 2 - 6c2 . (A54) 



These solutions in fact apply for essentially all gauges 
until the epoch of dark energy domination. The only 
exceptions are those that place explicit conditions on the 
dark energy fluctuations such as 5q — or uq = 0. 

It is instructive to consider the two limiting cases of 
Hubble drag domination and negligible Hubble drag or 
kinetic energy domination in Eqn. IJA36J1 . In the former 
case, 



p = 0, 



1 + 3wt 



(-2 + c\ + clw T ) • 



(A55) 



In terms of the field variables, the constant mode cor- 
responds to an initial condition where SX = 0, i.e. the 
kinetic energy terms in Eqn. (|A24p for the density per- 
turbation vanish. Therefore the potential energy fluctu- 
ation remains constant. Note that in spite of this the 
momentum density does not vanish. The second mode 
corresponds to an initial conditions with comparable ki- 
netic and potential energy (SQ terms) in the perturba- 
tion. An arbitrary initial condition can be decomposed 
into a superposition of the modes. Note that the second 
mode is decaying for c\ < 2/(1 + Wt) and growing for 
Cg > 2/(1+wt)- For a canonical kinetic term, this is a de- 
caying mode for all wt < 1- Nonetheless, if 2 > c\ > 3/2 
a small initial density or field fluctuation will be am- 
plified during radiation domination and freeze in during 
matter domination even though the background field is 
potential energy dominated during the whole expansion 
history wq rs — 1. 

The existence of this mode could potentially allow a 
solution to the gravitational wave problem of §V\ by am- 
plifying the field fluctuations from their initial conditions. 
However the amount of amplification is dependent on the 
initial ratio of kinetic to potential energy in the fluctua- 
tion. It also depends on the ratio in the background since 
the fractional fluctuation in the kinetic energy density (as 
opposed to the total energy density) must also remain 
small for the mode analysis to remain valid. Finding an 
explicit model that satisfies these conditions is beyond 
the scope of this work. 

In the opposite regime of kinetic energy domination, 
the two solutions become 



p = 0, 



3{2c 2 e - 1 - w T ) 

1 + 3wt 



(A56) 



The p — solution corresponds to a kinetic energy dom- 
inated perturbation where SX scales with X. The other 
solution corresponds to a pure velocity isocurvature mode 
where UQ/{krj) > iq ~ 0. Note that uq grows if 
Cg > 1/3 as is the case for the canonical kinetic term. 
In the field representation, this mode represents a case 
where the density perturbation is dominated by the po- 
tential energy and hence negligible in the kinetic energy 
dominated regime. The field fluctuation SQ then remains 
constant but the momentum density uq in Eqn. (|A27|) 
grows due to the redshifting of pq in the background. If 
the field later exits from kinetic energy domination, the 
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field fluctuation then becomes an energy density fluctu- 
ation. 

Finally, there is the well-studied tracking regime where 
1 + wq oc 1 + wt ■ Then 1 + wq is approximately constant 
and Eqn. I|A13|I gives (? a sw wq. The two solutions become 



2(1 + 3w T ) 



( 1 — 2wq + wt 



± 



^/(l + 2w Q + w T y - 8c2(l + w T )) ■ (A57) 



The index p is maximized by maximizing wq. The 
maximum wq — (? e for a kinetic energy dominated 
field and hence the fastest growing modes are given by 
Eqn. I|A56(1 . By further requiring the tracker condition 
wq < (1 + wt)/2 J43, we see that Re(p) < and there 



are again no growing modes. The field fluctuations then 
"track" and lose their dependence on the initial isocur- 
vature perturbation as is well known. Note that for the 
canonical kinetic term in the tracking regime, the pro- 
portionality is |45j 



l+WQ 



1 + wt 
2T-1 



(A58) 



where T was defined in Eqn. (|ll|l . The perfect tracker is 
attained at the limiting case of 1 + wq = 1 + wt or T = 
1 which is acheived for an purely exponential potential 
Vq cx e~ c< 2 where C is constant |63|. Here the dark 
energy remains a constant fraction of the total energy 
density. 
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